function main() {
	let scale=255.0/5000 ;
	let ds=pe.Dataset('modis/mod09cmg', 20230901000000);
	let ds1 = ds.forEachPixel( function(v,i) {
		r=v[0];//modis band1 for red
		g=v[3];//modis band4 for green
		b=v[2];//modis band3 for blue
		a=0;
		if( r>0 ) {
			r=r*scale;
			g=g*scale;
			b=b*scale;
			a=255;
		}else{
			r=g=b=0;
		}
		return [r,g,b,a] ;
	}) ;
	let ds2=ds1.clip3('canada.geojson',0);
	return ds2;
}